Thermal and electromagnetic properties of 166 Er and lfci7 Er 



167i 



E. Melby, M. Guttormsen, J. Rekstad, A. Schiller, and S. Siem 
Department of Physics, University of Oslo, N-0316 Oslo, Norway 



o 
o 
o 

(N 

> 

o 
O 
o 

CO 



> 

On 

o 
o 

o 
o 

">< 

<D 
i 

=5 



x 



A. Voinov 

Frank Laboratory of Neutron Physics, Joint Institute of Nuclear Research, 141980 Dubna, Moscow reg., Russia 

The primary 7-ray spectra of 166 Er and 167 Er are deduced from the ( 3 He,Q7) and ( 3 He, 3 He'7) 
reaction, respectively, enabling a simultaneous extraction of the level density and the 7-ray strength 
function. Entropy, temperature and heat capacity are deduced from the level density within the 
micro-canonical and the canonical ensemble, displaying signals of a phase-like transition from the 
pair-correlated ground state to an uncorrelated state at T c ~ 0.5 MeV. The 7-ray strength function 
displays a bump around _E 7 ~ 3 MeV, interpreted as the pygmy resonance. 

PACS number(s): 21.10-Ma, 24.10.Pa, 24.30. Gd, 27.70. +q 



I. INTRODUCTION 



The energy distribution of primary 7-rays provides 
information on both the level density and the 7-ray 
strength function. The present aim is to study the level 
density and 7-strength function of the nuclei 166 > 167 Er 
and to compare with other nuclei in this mass region. In 
addition, thermodynamic and electromagnetic properties 
of the two nuclei will be investigated. 

Most of the experimental information on the level den- 
sity of rare earth nuclei originates from the neutron- 
resonance spacing at the neutron-separation energy |l]] 
and direct counting of levels in the vicinity of the ground 
state Q| . Experimental knowledge of the 7-strength func- 
tion is mainly based on the study of photo-absorption 
cross sections ||] and radiative neutron capture Q] . 

A new method derives the level density and 7- 
strength function simultaneously from a set of primary 
7 spectra. This technique has proven to give a valuable 
supplement to the more traditional methods for level den- 
sity extraction, and to our knowledge, it represents the 
least model-dependent method to obtain the 7 strength 
function over a wide energy region below the neutron- 
separation energy. 

The nuclear level density is closely related to the ther- 
modynamic properties of nuclei, and can therefore be uti- 
lized to deduce e.g. entropy, temperature and heat capac- 
ity of nuclei. The thermodynamic observables may reveal 
phase transitions. 

The 7-strength function is a measure for the average 
electromagnetic properties of nuclei and has a fundamen- 
tal importance for the understanding of nuclear structure 
and reactions involving 7 rays At ~ 3 MeV of 7 
energy, a bump is observed in the 7-strength function 
of rare earth nuclei from ( 3 He,a) experiments [Q, and 
is probably of the same origin as the pygmy resonance 
found in (71,7) reactions [||. 

Section II describes the experimental methods. In 
Sect. Ill the experimental level density and 7-strength 
function of 166 Er and 167 Er are obtained. Section IV 



examines thermodynamic properties within the micro- 
canonical and the canonical ensemble, while electromag- 
netic properties of the two nuclei are discussed in Sect. V. 
Concluding remarks are given in Sect. VI. 



II. EXPERIMENTAL METHOD AND 
TECHNIQUES 

The experiment was carried out at the Oslo Cyclotron 
Laboratory. The reactions employed were the neutron 
pick-up ( 3 He,ev7) and the inelastic scattering ( 3 He, 3 He'7) 
with a beam energy of 45 MeV, populating 166 Er and 

Er with high excitation energy and low spin. The 
Q- value of the ( 3 He,a)-reaction is 14.142 MeV. The 
self-supporting 167 Er target was isotopically enriched to 
95.6% and had a thickness of 1.5 mg/cm 2 . 

The charged particles and 7 rays were recorded with 
the detector array CACTUS || , which contains 8 Si par- 
ticle telescopes and 27 Nal 7-ray detectors. Each tele- 
scope is placed at an angle of 45° relative to the beam 
axis, and comprises one Si front and one Si(Li) back de- 
tector with thickness 140 and 3000 /im, respectively. The 
Nal 7-detector array, having a resolution of ~ 6 % at 7 
energy E 1 — 1 MeV and a total efficiency of ~ 15 %, sur- 
rounds the target and particle detectors. In addition, two 
Ge detectors were used to monitor the spin distribution 
and selectivity of the reactions. 

The excitation energy of the resulting 166 Er and 167 Er 
nuclei are determined by means of reaction kinematics 
of the ejectile. By setting proper gates in the particle 
spectra, each coincident 7 ray can be assigned to a 7 cas- 
cade originating from a specific excitation energy. The 
data are sorted into a matrix of (E,E^) energy pairs, E 
being the excitation energy. Examples of the recorded 7 
spectra, the so called raw spectra, from two excitation 
energies are shown in the left panel of Figs. |l| and || for 
166 Er and 167 Er, respectively. Note that the statistics 
of 167 Er is about twice as good as for 166 Er, since the 
( 3 He, 3 He') reaction has a higher cross section than the 
( 3 He,a) reaction. 
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In order to determine the true 7-energy distribution, 
the 7 spectra are corrected for the response of the Nal 
detectors with the unfolding procedure of Ref. [|o| . Un- 
folded 7 spectra are shown in the central panel of Figs, [j] 
and H for the two nuclei. 

The now corrected (E,E 7 ) matrix comprises the 7- 
energy distribution of the total 7 cascade, and make it 
possible to derive the primary 7 matrix according to the 
subtraction technique of Refs. pl]Jl^ |. Primary 7 spectra 
can be seen in the right panel of Figs. [I] and g. 

The method of extracting the primary 7 spectra is 
based on the assumption that the decay properties of 
a bin of excited states are independent of whether the 
states are directly populated through the nuclear reaction 
or from de-excitation from higher excited states. This is 
believed to be approximately fulfilled because the long 
life time of excited states give the nucleus time to ther- 
malize prior to the 7 decay. 

III. EXTRACTION OF LEVEL DENSITY AND 
7- RAY STRENGTH FUNCTION 

The primary 7 matrix provides information on both 
the level density and the 7-ray strength function, en- 
abling a simultaneous determination of the two functions. 
The fundamental assumption behind the extraction pro- 
cedure is the Brink- Axel hypothesis ]l3| , p^ |, where the 
probability of 7-decay in the statistical regime, repre- 
sented by the primary 7 matrix P(Ei,E~), can be ex- 
pressed simply as a product of the final-state level density 
p(E — £L) and a 7-energy dependent factor F{E 1 ) 

P{E,E 1 )(xF(E 1 )p{E-E 1 ). (1) 

The 7-energy dependent factor is proportional to 

X; X L^7 i+1 /xL(£ 7 ), wner e Ixl(E 7 ) is the 7-ray 
strength function for the multipolarity XL. 

During the last few years a method for simultaneous 
deduction of level density and 7-strength function from 
primary 7 spectra has been developed at the Oslo Cy- 
clotron Laborato ry p |, ^2p5| ] . In the previously published 
results on 166 Er Jig]" the level densities and 7-strength 
functions were still unnormalized. The procedure is now 
fully replaced by an iteration technique deriving the level 
density and the 7-energy dependent factor by a x 2_n t to 
the experimental primary 7 spectra [El . This technique is 
described in detail in Ref. j^] , and only an outline of the 
normalization procedures is given here. It is important to 
note that the extracted experimental quantity p is in fact 
the density of levels accessible to the nuclear system in 
the 7-decay process. This quantity is interpreted as the 
level density, but may be influenced by selection rules in 
the 7 decay. For a thermalized nucleus in the continuum, 
the interpretation is approximately valid jlTj . 

Equation ([j]) has an infinite number of solutions. It can 
be shown || that all equally good solutions of Eq. (Q) can 
be obtained by the transformations 



p(E - E y ) = Aexp[a(E - E 7 )]p{E - E y ), (2) 
FiE^) = Bexp(a J B 7 )F(£; 7 ) (3) 

of any particular solution (p, F) . Consequently, neither 
the slope nor the absolute value of the two functions 
can be obtained through the iteration procedure, but the 
three variables A, B and a of Eqs. (g) and (||) have to 
be determined independently to give the best physical 
solution of the level density and 7-strength function. 



A. The level density 

The parameters A and a of Eq. (|^) can be determined 
by fitting the level density from the iteration procedure to 
the number of known discrete levels at low excitation 
energy and to the level density estimated from neutron- 
resonance spacing data at high excitation energy Q . This 
normalization procedure is shown for 167 Er in the right 
panel of Fig. ||. The deduced level density is fitted to the 
discrete levels (between the arrows in the upper panel of 
Fig. [3]) as far up in energy as we can assume that all levels 
arc known. At high excitation energies the deduced level 
density is fitted (between the arrows in the lower panel of 
Fig. ||) to a Fermi-gas approximation of the level density 
(line) forced to pass through the level-density estimate 
at the neutron-separation energy (filled square) obtained 
from the neutron-resonance spacing data. 

Unfortunately, 165 Er is an unstable nucleus (Ti/2 = 
10.36 h), making the neutron capture 165 Er(n,7) 166 Er 
difficult to investigate. Thus, there is no neutron- 
resonance data available for 166 Er. In this case the level 
density at the neutron-separation energy is estimated 
with the Fermi-gas expression 

exp(2V^7) 
12\/2a 1 /4[/5/4 (T ' ^ > 

where the level density parameter a = 0.21A 087 , the 
shifted excitation energy U = E — C\ — A, and the spin- 
cutoff parameter a is defined by er 2 = 0.0888\/atL4 2 / 3 , 
while A is the mass number, A the pairing-gap parame- 
ter and C\ = — 6.6A~ - 32 . All parameters are calculated 
according to Ref. which provides a level density in 
acceptable agreement with the experimental level density 
at the neutron-separation energy for 167 Er. Also the nu- 
clei 161 ' 162 Dy and 171 > 172 Yb show conformity between the 
experimental and theoretical level density at the neutron- 
separation energy |l9| . The left panel of Fig. || shows the 
normalization of 16 ^r, where the open square represents 
the theoretically calculated level density at the neutron- 
separation energy with an error of 80% jl8| . 

The normalized level density of 166 Er is compared to 
the level densities of 162 Dy and 172 Yb in the left panel 
of Fig. |J. The three level-density functions coincide well 
at high excitation energies, indicating that the chosen 
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value for the level density at the neutron-separation en- 
ergy for 166 Er is reasonable. For excitation energies below 
1.9 MeV the experimental level density of the erbium nu- 
cleus is seen to be lower than for the two others. In this 
energy region most levels in the three nuclei are believed 
to be known. From counting of known levels it is 
found that between 1.1 and 1.9 MeV, the level density of 
162 Dy and 172 Yb indeed exceed the level density of 166 Er 
with the same degree as observed experimentally. 

The right panel of Fig. || compares the level density 
of 167 Er with the level densities of 161 Dy and 171 Yb @. 
The two latter nuclei were populated through ( 3 He,a) 
reactions. The 171 Yb-level density is seen to be lower 
than the two others below 1.1 MeV. According to 
this agrees with reality below 900 keV. Above this limit 
the level densities from counting of known levels start to 
decrease, making it evident that not all levels are known. 

In Ref. [jllj the level density of 162 Dy extracted from 
two different reaction mechanisms are compared. The 
primary 7 spectra from the ( 3 He,a) pick-up reaction and 
the inelastic (' ! Hc.' ! He') reaction are found to provide 
level densities with very good agreement except below 
~ 1 MeV of excitation energy. There, the level density 
obtained from the inelastic 3 He scattering overestimates 
the number of levels by a factor of ~ 3, because this re- 
action populates more collective excitations with similar 
structure as the ground-state band. This causes a rela- 
tively large 7-transition rate from the direct populated 
states to the ground band, which in turn gives an overes- 
timation of the level density there. The same is observed 
in the level density of 172 Yb deduced from the inelas- 
tic reaction )2Cj ]. This effect of overestimation can not 
be ruled out in the level density of 167 Er derived from 
the ( 3 He, 3 He') reaction, however, the experimental data 
fit the known level density at low excitation energy very 
well, see Fig. ||. A difference between these ( 3 He, 3 He') 
experiments is the spin of the target nuclei. The 167 Er 
target has a spin of 7/2 making excitations to a large 
number of states probable, with a resulting fragmenta- 
tion of collective strength and a lower probability for 
subsequent 7 transitions directly into the ground-state 
band than in the even-even nuclei. Thus, the accessi- 
ble states in the 7 decay seem to approximates the level 
density quite well. 



B. The 7-ray strength function 



Blatt and Weisskopf |21| suggested the ratio of the 
partial radiative width Ti(Ey) and the level spacing of 
the initial states i with equal spin and parity Di to de- 
scribe the 7 decay in the continuum. The correspond- 
ing definition of the 7-ray strength function is given by 
I xl = T i (E 1 )/(E^ L+1 D i ), where X denotes the electric 
or magnetic character, and L defines the multipolarity of 
the 7 transition. 

After the normalization of the level density, the pa- 



rameter a of Eq. (H) is already fixed, and the slope, i.e. 
the exponential exp(ai? 7 ), is included in the 7-energy de- 
pendent functions shown for 166 Er and 167 Er in Fig. ||. 
Still, the normalization constant B of Eq. (||) remains to 
be determined. 

The 7-energy dependent factor F{E 1 ) is proportional 
to the sum of E^ L+1 fxL(E 7 ) for all possibilities of X 
and L. We assume that the 7 decay in the continuum of 
nuclei with low spin is mainly governed by electric and 
magnetic dipole radiation and that the accessible levels 
have equal numbers of positive and negative parity states. 
Thus, the observed F can be expressed by a sum of the 
El and Ml 7-strength functions only: 

BF(E^ - [/ei(£ t ) + /mi(£ 7 )]£ 3 . (5) 

The average total radiative width of neutron resonances 
0^7 ) H with excitation energy equal to the neutron- 
separation energy S n , spin /, and parity ir 



<r 7 > = 



p[t> n ,i,n) XLIf 7rf Jo 

p{S n - E 7 ,I f ,ir f ), (6) 

can be written in terms of F by means of Eq. (§). With 
the experimental level density p already normalized, the 
normalization constant B can be deduced. Assuming s- 
neutron capture, J, tt take the values It ± ^ , 7r f , where 
It , 7r t are the spin and parity of the target nucleus in the 
(n, 7) reaction. The experimental value of the total ra- 
diative width is the average over the possible spins of the 
compound state (I t ± i). Since we expect the 7 decay 
to be governed by dipole transitions, the second sum is 
restricted to possible final states with spin If and parity 
-ir/ accessible by dipole radiation. 

The methodical difficulties in the primary 7-ray extrac- 
tion, prevents determination of the functions F{E 1 ) and 
p(E) in the interval E 1 < 1 MeV and E > S n - 1 MeV, 
respectively. In addition, the data at the highest 7- 
energies, above E~ t ~ S n — 1 MeV, suffer from poor statis- 
tics. Therefore, extrapolations of F and p were necessary 
in order to calculate the integral of Eq. (|^). For the level 
density a Fermi-gas extrapolation is used (see Fig. ||), 
and for the 7-energy dependent factor a pure exponen- 
tial of the form exp(M? 7 ) is utilized (see Fig. ||). The 
contribution of the extrapolation to the total radiative 
width in Eq. @ does not exceed 15%, thus the errors 
due to a possibly poor extrapolation are expected to be 
of minor importance pj . 

Values for the neutron-resonance radiative width are 
given in |§, and for 167 Er we find = 92(8) meV. 
Since the radiative width of 166 Er is unknown, its value is 
taken as the average of (T 7 ) of the neighbouring isotopes. 
This is justified by the general slow variance in the (F 7 ) 
values of neighbouring isotopes for other elements. From 
interpolation, we thus adopt the value (r 7 ) = 90(20) 
meV for 166 Er. The normalized 7-strength functions of 
166 Er and 167 Er are compared to the strength functions 
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of 162 Dy and 172 Yb, and 161 Dy and 171 Yb, respectively, 
in Fig. |6[ All strength functions are increasing smoothly 
with 7 energy except for a bump around E 1 ~ 3 MeV. 
The 7-strength functions of neighbouring isotopes display 
a striking qualitative similarity. This can be expected for 
nuclei with approximately the same charge distribution. 
The location of the bump can however be seen to move 
towards higher 7 energies with increasing mass number. 



IV. THERMODYNAMIC PROPERTIES 

For a statistical description of hot nuclei, the micro- 
canonical thermodynamics is the proper theory. Within 
this frame, the system is isolated, giving a well de- 
fined energy. However, the canonical ensemble, permit- 
ting heat exchange, and the grand-canonical ensemble, 
which in addition allows particle exchange, arc often used 
due to mathematical difficulties with detailed calcula- 
tions within the micro-canonical ensemble. In this work, 
both the micro-canonical and the canonical ensemble will 
be utilized to discuss thermodynamic properties of the 
166,167 Er nuclei experimentally. 

In particular we will focus on different ways, within 
the two ensembles, to obtain experimental values of the 
critical temperature for the pair-breaking process and 
the general quenching of pair correlations. The micro- 
canonical ensemble gives a detailed description of the 
breaking of one, two, three,... nucleon pairs as a function 
of excitation energy, while the canonical ensemble reveals 
the general average properties of this phase-like tran- 
sition. The proper definition of thermally driven first- 
and second-order phase transitions in systems with few 
particles is a long standing problem, which will not be 
discussed in the present experimental work. 



A. Micro-canonical ensemble 

The micro-canonical partition function is simply the 
multiplicity of nuclear states, which experimentally cor- 
responds to the level density of accessible states. Thus, 
the experimental level density p(E) is our starting point 
for the extraction of thermodynamic properties of nuclei. 

The entropy is determined by 

S(E) = So + In P (E), (7) 

where the Boltzmann constant for simplicity is set to 
unity (ks = 1), and the normalization constant So can 
be adjusted to fulfill the condition of the third law of ther- 
modynamics; S — > when T — > 0, T being the nuclear 
temperature 



The ground-state band of even-even nuclei is assumed to 
have T = 0. Therefore, So is determined so that the en- 
tropy of the ground-state band in 166 Er is approximately 
zero. 

Figures ^ and || show the entropy and the temperature, 
respectively, deduced in the micro-canonical ensemble for 
166 Er and 167 Er. The entropy curve plotted in a linear 
scale is essentially identical to the level-density curve in 
a logarithmic scale. The small bumps in the entropy 
curves of Fig. ^ are enhanced through the differentiation 
performed in Eq. (||) to obtain the temperatures of Fig. ||. 

The specific heat can be determined from differentiat- 
ing the temperature 

<w>-(sg)~' (9) 

and is displayed in the lower panel of Fig. || for 166 Er 
and 167 Er. The double differentiation of the entropy 
has introduced strong fluctuations in the specific heat. 
Still it is distinct that the oscillations in the temperature 
give rise to successive positive and negative heat capac- 
ity. The spectacular feature of negative heat capacity is a 
direct consequence of the decrease in the micro-canonical 
temperature, and has recently been observed experimen- 
tally in the critical region of nuclear fragmentation in Au 
quasi-projectile sources formed in Au+Au collisions p3| . 

An anomalous decrease in the temperature means an 
unusual increase in the entropy, i.e. a more than normal 
opening of new domains of the phase space, e.g. open- 
ing of new degrees of freedom as more particlepairs are 
broken. The successive bumps of Figs. [7] and || indicate 
that the transition from a paired to an unpaired phase is 
a gradual process of breaking up more and more nucleon 
pairs. 

When we interpret the area of negative heat capacity, 
corresponding to the region where the temperature has 
a negative slope, as the location of the break up of a nu- 
cleon pair, a signal of the break up of the first pair can be 
seen around E ~ 1.5 MeV in 167 Er in Fig. |[ This value 
is close to twice the pairing-gap parameters 2A„ = 1.7 
MeV and 2A p = 1.5 MeV, which are the expected cost 
of breaking of a neutron and proton pair, respectively. 
The pairing-gap parameters are calculated from the em- 
pirical masses of a sequence of isotopes or isotones Q 
of 167 Er. The first negative slope in the temperature as 
a function of excitation energy in 161 Dy, 162 Dy, 171 Yb, 
and 172 Yb J25| also coincide roughly with 2A. A similar 
argument is however less successful for 166 Er. Deviations 
from 2A in the localization of the break up of the first 
pair of particles can be due to the influence of several 
structural effects in the nuclei, such as e.g. the Fermi- 
level position in the Nilsson single-particle scheme, vari- 
ation in the density of single-particle orbitals, and two 
quasi-particle couplings to collective degrees of freedom. 

It has been shown |2(| that phase transitions of finite- 
size systems can be observed in the micro-canonical en- 
semble without invoking the thermodynamic limit. Ac- 
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cording to |27|], a phase transition of first order can be 
recognized from the caloric curve in the region of co- 
existing phases, i.e. in the excitation energy region where 
the phase transition takes place. There, both the phase 
the system is leaving and the phase the system is attend- 
ing are existing simultaneously. For this co-cxistcncy the 
system has to "pay" an amount of entropy which is even- 
tually returned as the whole system is converted into the 
new phase Q. The non-oriented area between the in- 
verse temperature T(E)^ 1 and the line T c _1 is twice this 
entropy. The critical temperature has to satisfy the con- 
dition that the oriented area between T(E)^ 1 and T c _1 
equals zero. 

The process studied in the present work, namely the 
breaking of nucleon pairs and the general quenching of 
pair correlations, deviates from the process of multifrag- 
mentation |26| , p7|| . In the latter process, the geometrical 
surfaces of the fragments play a crucial role. In the de- 
pairing process, the long-range two- nucleon force has cor- 
relation lengths longer than the diameter of the nucleus. 
Thus, the phase-like situation of having a pair coupled 
or not, reveals no geometrical surface of interaction. 

In order to qualitatively analyze our experimental find- 
ings along the procedure of Ref. J2?J , the excitation en- 
ergy region around the outstanding bump in the temper- 
ature of 167 Er is expanded, and the entropy, inverse tem- 
perature and heat capacity of 167 Er are plotted in FigjE . 
The qualitative similarity with the picture of Ref. |27 ] 
suggests that the general process of breaking a nucleon 
pair around E ~ 2A appears like a phase transition of 
first order. 

Accepting this theoretical result, however, the similar 
experimental signals of breaking of further nucleon pairs 
indicate that the process of quenching of pairing corre- 
lations is a series of first order phase transitions. This 
interpretation might seem physically unattractive, and 
therefore we are anticipating further theoretical clarifica- 
tions. 

The convex shape of the entropy curve of Fig. ^ in- 
dicates a transition from the one quasi-particle regime 
to the three quasi-particle regime. The transition takes 
place between the two tangent points of the entropy and 
its enveloping curve (dashed line in the entropy plot). 
The transition region is marked with two dotted vertical 
lines. The critical temperature for breaking of the first 
nucleon pair in 167 Er is derived to be T c — 0.51(4) MeV 
(horizontal dashed line in the caloric curve of Fig. ||). 

The inter-phase entropy 1 could in principle be ex- 
pected to be equal to the envelope entropy, that is, an 
interpolation of the entropies before and after the pair- 
breaking process. But compared to this expectation, 
the entropy is seen to be lowered by ~ 0.6 units, cor- 
responding to the sum of the primary loss and the final 



1 As stated above, these phases are not physically separated. 



gain in entropy through the transition. It is unclear how 
this entropy-mixing contribution shall be interpreted in 
a small nuclear system. 

The first law of thermodynamics states that supplied 
heat equals the change in internal energy plus the work 
performed by the system: dQ = dE + dW. Since the 
nucleus is an isolated system the work dW is zero, and 
the increase in energy equals the supplied heat through 
the phase transition. In macroscopic systems, where the 
temperature is constant through a phase transition of 
first order, this increase in energy is the latent heat. We 
interpret the transition-region energy of approximately 
1.8 MeV (see Fig. |^) as the latent heat, even though the 
temperature is not constant. We note that this value 
is also very close to the expected cost 2A, of breaking 
up a particle pair. The dashed lines in the specific heat 
are purely meant as guiding lines towards the theoretical 
prediction of poles between the successive positive and 
negative branches. 

The derivative of the entropy-envelope curve in Fig. || 
equals the critical temperature. The somewhat stronger 
concave shape of the micro-canonical entropy of 167 Er 
compared to 166 Er (see Fig. [?]) suggests that increasing 
temperature is necessary to break further nucleon pairs 
of 167 Er, whereas the critical temperature for breaking of 
particle pairs is rather constant in 166 Er. After the break 
up of the first particle pair in 167 Er, the slope of the en- 
tropy becomes however more constant. A linear fit to the 
entropy above E ~ 2 MeV gives a critical temperature 
of 0.53(10) MeV for 166 Er and 0.55(4) MeV for 167 Er, 
indicating that the critical temperature for breaking of 
nucleon pairs in 167 Er has a slight increase. This can 
be understood from the blocking effect of the unpaired 
nucleon, increasing the distance to the Fermi surface for 
low-lying orbitals with coupled pairs. 

The critical temperature for breaking of the first par- 
ticle pair in 167 Er could be nicely determined from equal 
areas in the caloric plot. Unfortunately, the thermody- 
namic quantities of small systems are strongly fluctuat- 
ing. Thus, the critical temperature for each back bend 
of the temperature, e.i. breaking of a particle pair, can 
generally be difficult to obtain by this method. The slope 
of the entropy-envelope curve can however always serve 
as a measure. 



B. Canonical ensemble 

The transformation from the micro-canonical to the 
canonical ensemble is performed by the canonical parti- 
tion function 

oo 

Z(T) = J2"(E)e~ E/T - (10) 

E=0 
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The partition function is a Laplace transform of the mul- 
tiplicity of states lo(E) — AEp(E), where p(E) is the 
level density of accessible states in the nuclear reaction 
at the discrete energy E given in energy bins AE. 

The thermal average of the excitation energy in the 
canonical ensemble is 



(E(T)) = Z- 1 J2 Eco(E)e~ E/T . 



(11) 



E=0 



By the Laplace transform in Eq. ( J10| ) much of the infor- 
mation contained in the micro-canonical level density be- 
comes smeared out. This smoothing effect can be quan- 
tified by the standard deviation for the thermal average 
of the energy 



V(E 2 ) - {E)\ 



(12) 



which e.g. gives <j e = 3 MeV at E = 7 MeV. Thus, fine 
structure in the thermodynamic observables in the micro- 
canonical ensemble will not be visible in the canonical 
ensemble. The lines in the upper panel of Fig. || display 
the smooth variance of the canonical temperature as a 
function of the thermal average of the excitation energy. 

The mathematical justification of Eqs. (|l^) and Jll] ) 
is that the summation is performed to infinity. However, 
for the typical temperatures studied here, T < 1 MeV, 
the level density should be known up to 40 MeV, which 
was found [ ^8| to be a sufficiently high upper limit for 
the summation. The experimental level density is only 
covering the excitation energy region from zero to about 
S„ — 1 MeV. In the region of and above the neutron- 
separation energy S n , the Fermi-gas model is believed to 
describe the nuclear properties. Therefore, the experi- 
mental level density is extrapolated to higher excitation 
energies with the Fermi-gas approximation of Ref. [ p"8[ , 
see the solid lines of Fig. 

The heat capacity in the canonical ensemble is the 
derivative of the thermal average of the excitation en- 
ergy 



even 



d<E> 
dT ' 



(13) 



and it is shown for 166 . 167 Er as a function of the thermal 
average of the excitation energy in Fig. || and as a func- 
tion of temperature in Fig. [1^. The averaging done by 
the partition function can be seen to provide a smooth 
energy and temperature dependence of the heat capacity. 
The heat capacity of both nuclei shows an S shape as a 
function of temperature. This feature is interpreted f2§| l 
as a fingerprint of a phase transition from a state with 
strong pairing interaction to a state where the pairing 
correlations are quenched. Because of the smoothing 
performed by the canonical-partition function, discrete 
transitions between the different quasi-particle regimes, 
as observed within the micro-canonical ensemble, arc hid- 
den, and only the phase transition related to the quench- 
ing of the pair correlations as a whole can be seen. 



The shape of the heat-capacity curve is related to the 
level-density. In general a constant-temperature level 
density gives a pronounced S shape while a pure 
Fermi-gas level density provides a linear heat capacity 
C v = 2aT - 3/2; see the dashed lines of Fig. [l0| The 
more prominent S shape in 166 Er is due to a quicker rise 
in the heat capacity. This is the same phenomenon as 
observed in the micro-canonical ensemble, where the con- 
cave shape of the 167 Er entropy was assigned to the block- 
ing effect of the odd neutron. Shell model Monte Carlo 
calculations |2£j on various Fe isotopes have shown that 
the pairing-phase transition is strongly correlated with 
the suppression of neutron pairs with increasing temper- 
ature. It is also observed that the reduction of neutron 
pairs is much stronger in the even- than in the odd-mass 
isotopes, giving rise to a more pronounced S shape in the 
even nuclei. The same difference in the heat capacity is 
also observed experimentally between 161 Dy and 162 Dy, 
and 171 Yb and 172 Yb @. 

In order to extract a critical temperature for the 
quenching of pair correlations from the canonical data, 
two different approaches has been tested. First, the lo- 
cal maximum of the heat capacity relative to the Fermi- 
gas approximation of Fig. [l 7 ^ have been determined to 
be 0.47(10) MeV for 166 Er and 0.59(10) MeV for 167 Er, 
see the arrows of Fig. [Io| This method relies very 
much on the Fermi-gas extrapolation of the level den- 
sity. The second approach depends on the assumption 
that the nuclear level density can be approximated by 
a constant-temperature expression at low excitation en- 
ergies. Within the canonical ensemble this assumption 
provides the relations 



(E(T)Y 



C V (T) = (1 - T/t)- 



(14) 



(15) 



where r is the constant-temperature parameter. By fit- 
ting T -1 as a function of (E^))^ 1 to a straight line 
of slope 1, r can be determined. To minimize the de- 
pendency of the extracted critical temperature on the 
extrapolated level density, the fit is performed between 
(E(T)) - 0.5 - 2 MeV. According to Eq. (|lTj) this cor- 
responds to energies in the level density curve up to ~ 6 
MeV. Then r is identified with the critical temperature, 
since CV(T) has a pole at T = r. The critical tempera- 
ture is determined to T c = 0.44(10) MeV for 166 Er and 
T c = 0.52(9) MeV for 167 Er. The dashed-dotted lines of 
Fig. [n] describe Eq. ( |l5| ) with asymptotes at r = T c . 

The different values derived for T c are quoted in Ta- 
ble |. The micro-canonical caloric-curve method (Fig. ||) 
is only supposed to provide the critical temperature 
for breaking of the first particle pair, while the others 
should give the critical temperature for the total quench- 
ing of the pairing correlations. The canonical constant- 
temperature approach relies on experimental data, but 
it can be questionable whether Eqs. ([l4]) and ([if]) de- 
scribes the nucleus well at low excitation energies. The 



G 



method of deducing T c from the local maximum of the 
canonical heat capacity relative to the Fermi-gas approx- 
imation depends highly on the extrapolation of the level 
density. The most straight forward method is to obtain 
the temperature directly from the slope of the micro- 
canonical entropy. There seem to be a tendency for the 
critical temperature in 166 Er to be slightly lower than in 
167 Er. All results are however equal within the approx- 
imated uncertainties, with a mean value of 0.5(1) MeV 
for both nuclei. 

The results for 166 Er agrees with recent finite- 
temperature Hartree Fock Bogoliubov calculations for 
164 Er |30|], in which the critical-temperature estimate ex- 
tracted is 0.7 MeV. This value is somewhat higher than 
the present derived result for 166 Er. This diversity is 
probably due to application of different interpretations 
of where exactly the phase transition takes place. Since 
the shapes of the heat-capacity curves of 164 Er and 166 Er 
are practically identical, we assume that the calculated 
critical temperature for quenching of pair correlations in 
16 Er coincides with the experimental findings for 166 Er. 

Relativistic Hartree-BCS calculations |3lJ also find 
the proton and neutron pairing gaps to vanish around 
0.4 - 0.5 MeV of temperature in 166 Er and 170 Er. The 
heat capacity shows two peaks corresponding to these 
events, characteristic of second order phase transitions 
from superfluid to normal phase |31|. Thus theoretically 
the pair quenching is calculated to be a phase transition 
of second order. Provided that negative heat capacity is 
a sufficient signal, the observed breaking of nucleon pairs 
in the micro-canonical ensemble appears like first-order 
phase transitions. The process of breaking one nucleon 
pair does however not lead to a total quenching of the 
pair correlations, and may not be described by the same 
physics as the pairing-phase transition. It is not clear 
how a series of phase transitions should be interpreted 
physically. 



V. ELECTROMAGNETIC PROPERTIES 

The 7-strength function is governed by different multi- 
polarities of electric and magnetic character. We assume 
however, that the 7 decay in the continuum is dominated 
by dipole transitions and will try to model the 7-strength 
function theoretically. It is commonly adopted that the 
El strength is determined by the giant electric dipole 
resonance (GEDR) at high 7 energies. More doubtful is 
the assumption that the tail of the GEDR governs the 
El strength at low -renergies when the tail approaches 
zero @. A model (3^] taking into account the energy 
and temperature dependence of the GEDR width is of- 
ten utilized to describe the experimental data ||. The 
Ml strength function also plays an important role gov- 
erning the 7 emission for low 7 energies. Experiments 
indicate the existence of an Ml giant resonance due to 
spin- flip excitations in the nucleus ]33j]. 



In Fig. O, the experimental 7-strength function is fit- 
ted by a theoretical strength function taking into account 
both the giant electric dipole resonance and the spin-flip 
resonance. In addition, a weaker resonance at lower en- 
ergies is needed in order to fit the experimental data. Be- 
cause of the much lower strength compared to the GEDR, 
the resonance is denoted the pygmy resonance. To ac- 
count for the El radiation, the model fl32i 



Zei(£ 7 ) 



1 OJaEirl^+^T 2 ) 
(16) 



3ir 2 h 2 c 2 Eei(E% — E m ) 



is adopted. The values for the giant electric dipole reso- 
nance parameters ctei, Tei and E^i are taken from 0. 
The temperature parameter T is utilized as a constant 
fit parameter. The Ml radiation is described by 



cmi-E^T 2 ^ 



(17) 



where ctmi, Tmi and Emi are the giant magnetic 
dipole resonance parameters, which are taken from p2[ . 
The pygmy resonance is here described with a similar 
Lorentzian function / py as Eq. (jlT]), where the pygmy- 
resonance strength er py , width r py and centroid E py 
have been fitted in order to adjust the total theoretical 
strength function 



/ = K(f m + /j 



Ml, 



(18) 



to the experimental data. The resulting theoretical 7- 
strength functions are shown as solid lines in Fig. |ll|. The 
obtained values of the fitting parameters for the pygmy 
resonance, the normalization constant K and the tem- 
perature T are given in Table [0]. 

The low-energetic resonance has its centroid at E~ — 
2.98(8) MeV in 166 Er and at E 1 = 3.24(7) MeV in 167 Er. 
In (n, 7) reactions a pygmy resonance is observed in 
the energy region 2.5 — 3.5 MeV for nuclei of mass 160- 
176. Both the centroid and the width of the resonance 
is found Q to increase gradually with neutron number. 
In 166 Er the centroid, width and strength of the pygmy 
resonance are somewhat lower than in 167 Er, fitting very 
well into the systematics of @] . The enhanced 7 strength 
around 3 MeV in 161 . 162 Dy and 171 - 172 Yb § also fit into 
this systematics, indicating that the ~ 3 MeV bump af- 
ter ( 3 He,a) and ( 3 He, 3 He') reactions probably is of the 
same origin as the pygmy resonance reported in [|| . 

The pygmy resonance has been explained by the en- 
hancement of the El 7-strength function |||. The possi- 
bility that it on the contrary is of Ml character can not be 
excluded. At an excitation energy around 3 MeV, there 
is a concentration of orbital Ml strength in the weakly 
collective scissors mode [Q . The scissors mode was first 
observed in electron-scattering experiments pq] , and is 
confirmed by the (7,7') reaction ptj ]. A systematic sur- 
vey of centroids of observed Ml strength distributions in 
even-even nuclei in the A = 134— 196 mass region is given 
in Ref. |Bq|, where the scissors mode of 166 Er is found to 
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be centered at 2.961(26) MeV, in perfect agreement with 
the pygmy resonance found here. 

Figure |ll| also shows the predicted individual contri- 
butions from the giant electric dipole resonance /bi, the 
giant magnetic dipole resonance /mi and the pygmy reso- 
nance f py to the 7-strength function. The strength func- 
tion is generally dominated by El radiation. The Ml- 
strength function is always ~ 20% lower than the 
El-strength function fsi- Moreover, the pygmy reso- 
nance contributes considerably around 3 MeV. 

The normalization constant K is close to 1 for 167 Er 
(Table 0), showing that the adopted theoretical model 
reproduces the absolute values of the 7-strength func- 
tions very well. The deviation from 1 can be explained 
by uncertain values of the GEDR parameters for the in- 
vestigated nuclei and uncertainties in the experimental 
normalization of the 7-strength functions. Also the K 
value for 166 Er is reasonable, provided that the adopted 
value for the neutron-resonance radiative width is reli- 
able. 



VI. CONCLUSION 

Levels in 166 Er and 167 Er in the excitation region up 
to the neutron separation energy were populated with 
the ( 3 He,cry) and ( 3 He, 3 He'7) reaction, respectively. The 
level density and 7-ray strength function of 166 Er and 
167 Er are determined from their corresponding primary 
7-ray spectra. 

Thermodynamic observables are deduced from the 
level density and display signatures of phase-like transi- 
tions within the micro-canonical and the canonical en- 
semble, interpreted as the transition from a strongly 
pair-correlated phase to an uncorrelated phase. Micro- 
canonical thermodynamics give the possibility of investi- 
gating the successive breaking of nucleon pairs in detail, 
information which is hidden in the canonical approach. 
The canonical ensemble on the other hand, reveals the 
average properties of the pairing transition. In addition, 
the canonical ensemble provides an excellent opportunity 
to study the different mechanisms governing the ther- 
modynamic properties of odd and even systems. The 
increase in the heat capacity with temperature is much 
steeper in 166 Er than in 167 Er, probably due to the ten- 
dency of the odd neutron to block higher quasiparticlc 
excitations. Various estimates of the critical tempera- 
ture for the pairing-phase transition is performed, giving 
the result T c ~ 0.5 MeV. 

The experimental 7-strength function is fitted by a the- 
oretical strength function, assuming that the 7 decay in 
the continuum is governed by dipole transitions. The 
contribution of electric and magnetic dipole radiation to 
the 7-strength function is recognized. A bump is ob- 
served in the 7-strength function at 3.0 and 3.2 MeV in 
the 166 Er and 167 Er, respectively, and is probably of the 
same origin as the pygmy resonance found in the (n, 7) 



reaction. 
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TABLE I. Different extracted values for the critical temperature. 



Extraction method 


i66 Er 


167 Er 




T c (MeV) 


T c (MeV) 


micro-canonical caloric curve 




0.51(4) 


micro-canonical entropy 


0.53(10) 


0.55(4) 


canonical heat capacity 


0.47(10) 


0.59(10) 


canonical constant temperature 


0.44(10) 


0.51(9) 



TABLE II. Parameters obtained from the fitting of the 7-ray strength function. 



Nucleus 


E py 




1 py 


T 


K 




(MeV) 


(mb) 


(MeV) 


(MeV) 




ibb Er 


2.98(8) 


0.30(4) 


1.3(3) 


0.31(5) 


1.3(2) 


167 Er 


3.24(7) 


0.43(4) 


1.7(2) 


0.36(2) 


1.27(6) 
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FIG. 1. Gamma-ray spectra at E — 4 and 6 MeV of excitation energy for Er. Left panel: Raw 7-ray spectra, central 
panel: Unfolded 7-ray spectra, and right panel: Primary 7-ray spectra. 
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FIG. 2. Same as Fig. 1 for 167 Er. 
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FIG. 3. Normalization of the experimental level density (data points) of 166 Er (left panel) and 167 Er (right panel) between the 
arrows to known levels at low excitation energy (histograms) and to the Fermi-gas level density ( 166 Er) or to the level density 
calculated from neutron-resonance spacings ( 167 Er) at the neutron-separation energy (squares). The lines are the Fermi-gas 
approximations. 
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FIG. 4. The normalized level density of 166 Er (left) and 167 Er (right) compared to 162 Dy and 172 Yb, and 161 Dy and 171 Yb, 
respectively. The Dy isotopes are multiplied by 0.1 and the Yb isotopes by 10 for better visualization. All nuclei are popolated 
through the ( 3 He,a) reaction except for 167 Er, which is populated through the ( 3 He, 3 He') reaction. 
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FIG. 5. The 7-energy dependent factor of ' Er. The lines are extrapolations needed to calculate the normalization 
integral of Eq. (0) . 
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FIG. 6. The normalized 7-ray strength function of 166 Er (left) compared to 162 Dy and 172 Yb, and 167 Er (right) compared 
to 161 Dy and 171 Yb. The Dy isotopes are multiplied by 0.1 and the Yb isotopes by 10 for better visualization. 
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FIG. 7. The entropy as a function of excitation energy in Er and Er 
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FIG. 8. The temperature (upper panel) and heat capacity (lower panel) of 166 Er (left) and 167 Er (right) from the mi- 
cro-canonical ensemble (data points) and the canonical ensemble (lines). 
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FIG. 9. Entropy (upper panel), inverse temperature (central panel) and specific heat (lower panel) as functions of excitation 
energy in 167 Er. The vertical, dotted lines indicate the approximate region in which the phase-like transition takes place. See 
text. 
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FIG. 10. The semi-experimental heat capacity of 166 Er (left) and 167 Er (right) as a function of temperature. The dashed 
lines are the Fermi-gas approximations Cy = 2aT — 3/2. The local maxima relative to the Fermi-gas expressions are marked 
with arrows, while the dashed-dotted lines describe the estimates of Eq. (^|), where r is recognized as the critical temperature 
and marked with the vertival lines. 
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FIG. 11. The experimental 7-ray strength function (data points) of 166 Er (left) and 167 Er (right). The solid line is the fit to 
the data by the theoretical model. The dashed lines are the respective contributions of the GEDR, the GMDR, and the pygmy 
resonance to the total theoretical strength function. 
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